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Single molecule study of the DNA denaturation phase transition in the force-torsion 

space 
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We use the "magnetic tweezers" technique to reveal the structural transitions that DNA undergoes 
in the force-torsion space. In particular, we focus on regions corresponding to negative supercoiling. 
These regions are characterized by the formation of so-called denaturation bubbles, which have an 
essential role in the replication and transcription of DNA. We experimentally map the region of the 
force-torsion space where the denaturation takes place. We observe that large fluctuations in DNA 
extension occur at one of the boundaries of this region, i.e., when the formation of denaturation 
bubbles and of plectonemes are competing. To describe the experiments, we introduce a suitable 
extension of the classical model. The model correctly describes the position of the denaturation 
regions, the transition boundaries, and the measured values of the DNA extension fluctuations. 

PACS numbers: 82.37.Rs, 87.14.gk, 87.15.La 



The nanomechanics of DNA play an important role at 
the biological and biochemical levels jl|. Thus, under- 
standing the transcription and duplication phenomena is 
a relevant open topic to which a quantitative comprehen- 
sion of DNA mechanical characteristics is fundamental. 
In particular, because any transcription or duplication 
process implies the local and temporary separation of the 
two DNA strands (i.e., DNA breathing Q, y| or denat- 
uration bubbles 049), understanding denaturation rep- 
resents the first building block towards the theoretical 
comprehension of DNA metabolism. A well-known and 
promising technique for studying nanomechanical prop- 
erties is the magnetic tweezers (MT), which allows one 
to impose a stretching force and a torsion to a single 
DNA molecule while also monitoring the simultaneous 
extension of the same molecule [7|, |8[ . The versatility of 
the MT technique has been exploited to investigate DNA 
nanomechanics in the presence of proteins, enzymes, lig- 
ands, and drugs |0l— ll 31] , and phcnomcnologically analyzed 
[141] . The initial pioneering MT studies focused on the 
topology of DNA molecules and showed that torsion can 
produce a so-called "plectoneme" , which reduces DNA 
For modeling plectoneme formation, 



15, 1 



extension 

the DNA can be simply described as an elastic rod [17J . 
The experiments showed that the plectonemes disappear 
when the force becomes sufficiently high and the direc- 
tion of the torsion is toward the unwinding of the DNA 
double helix [8J. This chiral effect, which goes beyond 
the elastic rod model, has been_explained in terms of 
denaturation of the double helix 
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In this work we use the asymmetry between the DNA 
extension under positive and negative torsion as a hall- 
mark of denaturation. For the first time, we systemati- 
cally evaluate the occurrence of mechanical denaturation 
in the force-torsion space. We find that large temporal 
fluctuations of extension arise at one of the boundaries of 
the denaturation region. Finally, we interpret the experi- 
mental data with a simple mechanical model obtained by 



considering a denaturation term to the classical energy 



171 ] used to describe the buckling transition. 



Several MT apparatuses have already been reported 
in the literature and, in our set up [12|, we generally fol- 
low the most classically proposed schemes [3, |l9(. The 
technique is based on the following procedure: one end 
of the DNA is connected by standard biochemical tech- 
niques p to a commercial micron-sized superparamag- 
netic bead and the other DNA end is fixed to the inner 
wall of a squared capillary tube 2fJ, [2l( . Forces or tor- 
sions are then applied to the bead through a field gen- 
erated by external permanent mag nets, whose position 
and rotation can be controlled 22, 231. The movement 



of the magnetic bead is transferred to the DNA and thus, 
indirectly, forces or torsions are applied to the molecule 
The DNA extension L e is measured by considering 



24 



the diffraction images generated at different heights of 
the bead, which is illuminated by a LED light [19j. The 
force F exerted by the magnetic field on the bead and, as 
consequence exerted on the DNA, is obtained as in ref. 

Figure [1] shows the average value of the extension (L e ) 
of a DNA (~ 6000 base pairs) molecule as a function of 
the number of imposed turns n t at different values of the 
applied stretching force (F = 0.25 pN, 0.63 pN, 0.79 pN, 
and 1.14 pN). The data are the result of an averaging 
procedure on several values of L e (t) taken as a function 
of time at a frequency rate of 60 Hz during a time interval 
of several seconds. These results arc typical in the litera- 
ture of MT [l7j and are qualitatively described as follows. 
When increasing the absolute value of nt, the torsion is 
first absorbed by elastic twist deformations and the DNA 
extension remains approximately constant. Above the so 
called nb buckling transition, dependent on F and indi- 
cated by the vertical arrows in Fig. [1] the creation of 
plectonemes induces a progressive linear decrease in the 
DNA extension. At low forces (F < 0.6 pN), when in- 
creasing nt , the data show a symmetric trend at positive 




FIG. 1: Average DNA extension (L e ) as a function of the 
number of imposed turns nt- Data (dots) and theoretical 
curves (lines) are obtained for different values of the applied 
force F, as indicated in the figure. The short vertical arrows 
point to the buckling transitions n^. 
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FIG. 2: Asymmetry AL e (n t ) = (L e (-n t )) - (L e (n t )) (Panel 
A) and standard deviation of the DNA extension ctl,, (Panel 
B). measured as a function of the applied force F. Data ob- 
tained for different values of the imposed turns nt (indicated 
in the figures). Inset: measured maximum values aL e Max 
obtained from the <rz, e vs F curves shown in the main figure. 



and negative torsion. At high force values (F > 1 pN), 
the situation is different: at negative imposed turns, the 
L e vs nt curves are no longer symmetric, and the DNA 
extension is approximately constant due to the formation 
of so-called denaturation bubbles [18| . 

The asymmetric behavior of the (L e ) vs n t is quanti- 
fied by the expression AL e (n t ) = (L e (—n t )) — (L e (n t )). 
The AL e (n t ) values are shown in Fig. [2]A. as a function 
of F. We observe a transition between the plectonemic 
behavior at low forces and the formation of denaturation 
bubbles at high forces. This transition is highlighted by 
an increasing asymmetry and occurs around a character- 
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FIG. 3: Asymmetry AL e (n t ) = (L e (-n t )) - (L e (n t )) (Panel 

A) and standard deviation of the DNA extension gl c (Panel 

B) measured as a function of the imposed turns n t . Data 
obtained for different values of the applied force F (indicated 
in the figures). 



istic force F c h ar « 0.78 pN, which does not depend on 
n t . 

At the intermediate forces (0.6 pN < F < 0.9 pN) of 
Fig. [TJ we have also observed that the average extension 
(L e ) at negative values of n t (n t < -20) appears nois- 
ier than the data for a positive n t . We have quantified 
the fluctuations by calculating the standard deviation of 
the DNA extension <7Le — V '(Lf) — (L e ) 2 . We have also 
checked that ULe was not dependent on the total time of 
the average and that such time was much longer (more 
than 100 times) than the characteristic time of the corre- 
lation function of the data. The values of a^ e as a func- 
tion of F for various values of the number of imposed 
turns n t are shown in Fig. [5J3. An increase in the exten- 
sion fluctuations a^ e occurs in a narrow range of applied 



force and presents a maximum value ol, 



Ma 



around the 



characteristic force F rl 



the fluctuations arc associated 



with the denaturation transition. The value (JL e Max is 
an increasing function of \nt\. This is confirmed in the 
inset of Fig. [5J3, where the obtained values of o^^Max 
are plotted for different values of nt and show a linear 
behavior. 

We present the same quantities, AL e and ctlc as a 
function of n t for different values of force F (see Fig. [3J. 
We observe that the AL e increase with negative supcr- 
coiling starting at n t = -10 (panel A). This transition 
is not accompanied by an increase in o^e (panel B). In- 
stead, the strongest fluctuation takes place at the F c h a r, 
as shown previously. 

We will now present a simple mechanical model that 
can explain all the experimental findings, quantitatively 




FIG. 4: Panel A: Calculated phase diagram of DNA struc- 
ture in the (nt,F) space. White: DNA compatible with the 
WLC model. Blue: plectonemic regions. Green: zero DNA 
extension regions. Red: denaturation bubbles region. Yellow: 
coexistence region of denaturation bubbles and plectonemes 
at zero DNA extension. The vertical and horizontal lines 
correspond to the measurements presented at constant F (in 
Fig.[T]and Fig. [3} and at constant m (in Fig.[5J), respectively. 
Panel B (enlargement of panel A): contour plot of the cal- 
culated {level lines) olz- Crosses (circles) indicate measured 
values of <7Le < 0.15 /im (iri e > 0.15 /mi). 



confirming the existence of a fluctuation increase at a 
characteristic value of the force F c har and deriving an 
original relation between F c h ar and DNA nanomcchani- 
cal constants (bending constant B and the binding en- 
ergy between the DNA bases (denaturation energy)). 
The model is an extension of the classical theory of plec- 
toneme formation [171 ] considering the denaturation en- 
ergy. Indeed, when twisted by the magnetic bead, the 
DNA molecule relaxes the applied work in three dif- 
ferent ways: by twisting (twisting energy: Egoist), by 
forming plectonemes (plectonemic energy: E p i ect ) and 
by partially denaturating (denaturation energy: Ed en )- 
Accordingly [171 ]: 



1 C 

Efrwist = 77T _ (^ 7r ' n e) 



Epiect = (27Ti? • n p ) ■ (^2 + F ) 
while we have evaluated Eden as: 



Ede 



a ■ nd 



(1) 



(2) 



(3) 



number of plectonemes of radius R, rid is the number 
of turns relaxed by partial DNA denaturation, C is the 
twisting constant, and a is a phenomenological constant 
that will be discussed later. We impose the additional 
topological relation: n e = n t — n p — rid- Furthermore, we 
assume that the DNA extension variations are mainly 
due to the plectoneme formation: L e = Lq — 2-kR ■ n p . 
For this assumption, we disregard for the moment the 
Worm-Like-Chain (WLC) dependence of DNA extension 
on the applied force and we suppose that the denatu- 
ration phenomenon does not induce a significant DNA 
length variation. 

We can estimate the characteristic force Fehar by using 
the following back-of-the-envelope calculation. The char- 
acteristic force is obtained by equating the plectonemic 
energy and the denaturation energy which is necessary 
for relaxing to the same torsion angle (i.e. An p =Arid). 



The value R = -JB/2F of the plectonemic radius [17J is 
kept constant and obtained from the equilibrium of the 
twisting energy and the plectonemic energy, resulting in 
the following relations: 



(2nR) ■ ( 



B 

2~R2 



F) 



obtaining for the characteristic force: 



8Btt 2 



(4) 



(5) 



The parameter a is proportional to the binding energy 
between the bases of the two DNA single strands. Con- 
sidering the experimental value F c har ~ 0.78 pN and a 
DNA persistence length of about 50 nm, the resulting 
value of a is on the order of 10 -19 J, which is consistent 
with the reported average value of the free energy per 
DNA base pair, AG « 2.0 kcal/mole J25H27|. In a first 
approximation we consider the denaturation energy, de- 
scribed by the average single parameter a, to be indepen- 
dent on the specific DNA sequence. The resulting value 
of a corresponds to approximately 8 denaturated bases 
to relax one turn. This number of base pairs is consistent 
with the pitch of the double helix. More rigorously, fixing 
the external force F and the total imposed torsion n t , we 
can study the total energy landscape E to t as a function of 
the parameters n p , rid, and R. By minimizing the expres- 
sion of the total energy E to t in the (n t ,F) plane, we de- 
rive how DNA behaves under imposed torsion and force. 
The resulting calculations are presented in Fig.0]^ where 
the different colored regions represent different DNA be- 
haviors and the boundary lines correspond to transition 
lines. In the white central region, the DNA reaches the 
maximum extension compatible with the applied force 
according to the WLC model 17 1. Entering the lateral 



where Lq is the contour length of the DNA, n e is the 
number of turns that store energy in the twist, n p is the 



blue regions, plectonemic formation starts and DNA ex- 
tension is consequently reduced. The green regions corre- 
spond to achieving zero DNA extension. The novelty of 
the model is its capability for predicting the denaturation 



regions shown in the red (denaturation bubble) and yel- 
low (coexistence regions between plectonemes and bub- 
bles) rectangles. The boundary of each zone corresponds 
to a transition between different structural phases. We 
concentrate our attention on the transition to the denatu- 
ration state. It is possible to enter the denaturation zone 
(red and yellow) by crossing three different lines, but only 
the horizontal line delimiting the red zone is character- 
ized by a step in the DNA extension. Indeed, below this 
line the DNA has a reduced length due to the presence of 
several (depending on the turn position) plectonemes and 
above this line the DNA is fully extended. Accordingly, 
the fluctuation between the structures configuration is 
also accompanied by a oscillation in the DNA extension. 
This feature is unique to this transition; none of the other 
phase structure boundaries reported in Fig.HJ^ introduce 
any discontinuity in L e . 

Furthermore, because L e depends on n p , calculating 
n p using the proposed expression for the energy E to t and 
a Boltzmann distribution allows for the prediction of the 
average value of the DNA extension and its fluctuation 
(jLe- The calculated values of (L e ), corrected for the 
WLC model, are presented as a function of n t by the 
continuous lines in Fig. [I] showing a good agreement 
with the experimental data. The results are obtained 
assuming the following values: B = 1.29-10 -11 eV-cm, 
C = 3.62-ICT 11 eV-cm, L = 1-88 fjtm and a = 0.70 cV. 
As shown in Fig. [TJ the model accurately describes the 
classical plectonemic behavior observed for n t > 0, which 
was already well described in the past [17| . as well as 
the transition between the plectonemic and denaturation 
regime for n t < 0. 

As expected, the region characterized by non-negligible 
calculated DNA extension fluctuations 0£ e is located in 
a specific area of the plane (n t ,F): at the horizontal 
boundary of the red zone. This theoretically predicted 
region is shown as a contour plot in Fig. 3J3 (enlarge- 
ment of a region of Fig. UK). In Fig. 0J3 we also present, 
as crosses or circles, respectively, the measured standard 
deviations having negligible (o^e < 0.15 /im) or signif- 
icant (<7Le > 0.15 fim) values. From Fig. [4)3, we can 
appreciate the agreement between our model and the ex- 
perimental results: the regions of largest fluctuation are 
observed where the model predicts a marginal stability 
[28j | of L e , i.e., for low n t value (n t ^ -30) and near the 
characteristic value of the force F r j lar sa 0.78 pN. 

Moreover, in Fig. 2J\ we sketch vertical and horizontal 
thick black segments corresponding to the lines explored 
by the experiments presented in Fig. [5] (fixed n t ) and 
Fig. [3] (fixed F). As predicted by the model, the asym- 
metry in L e indicates entrance into the red zone and large 
fluctuations appear when crossing the red zone horizon- 
tal boundary and when the experiments are performed 
in its proximity. 

In conclusion, we have characterized the denaturation 
transition caused by external applied force and torsion, 



and we have introduced a nanomechanical model able to 
link the measured force at which the denaturation occurs 
to a parameter a related to the double strand stability. 
Because DNA denaturation and melting is at the origin of 
several biological problems ranging from DNA replication 
and transcription to the detection of transcription initi- 
ation points, the ability demonstrated here for a direct, 
quantitative single molecule measurement of the char- 
acteristic force and denaturation energy opens the way 
for several studies utilizing more sophisticated and bio- 
logically realistic situations and in the presence of DNA 
binding molecules or proteins [29J |. 
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